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Abstract. In this paper we present several examples of solving algorith¬ 
mic problems from the Google Code Jam programming contest with Pi- 
cat programming language using declarative techniques: constraint logic 
programming and tabled logic programming. In some cases the use of 
Picat simplifies the implementation compared to conventional impera¬ 
tive programming languages, while in others it allows to directly convert 
the problem statement into an efficiently solvable declarative problem 
specification without inventing an imperative algorithm. 


1 Introduction 

Google Code Jan0 (GCJ) is one of the biggest programming competitions in the 
world: almost 50,000 participants registered in 2014, and 25,462 of them solved 
at least one task. 

GCJ competitors can use any freely available programming language or sys¬ 
tem (including Picat 0 described in this paper). We show examples of solving 
GCJ problems with Picat using constraint logic programming and tabled logic 
programming. 

Picat is a new logic-based multi-paradigm programming language. Picat 
shares many features with Prolog, especially B-Prolog [4], but also has many 
distinct features: optional destructive assignments, functions in addition to pred¬ 
icates, explicit non-determinism, list comprehensions [5]. 

Picat implementation of TPK algorithm 

To give an idea of Picat syntax to a reader unfamiliar with the language we 
present an implementation of TPK algorithm. TPK is an algorithm proposed 
by D. E. Knuth and L. T. Pardo [2] used to show basic syntax of a programming 
language. The algorithm allows user to input 11 real numbers (ao ... aio). After 
that for i = 10... 0 (in that order) the algorithm computes value y = /(dj, 
where /(f) = i/|f[ + 5t 3 4 , and outputs a pair (i,y) if y < 400, or (*, TOO 
LARGE) otherwise. 

3 https : //code . google . com/code j am 

4 http://picat-lang.org/ 



1 f(T) = sqrt(abs(T)) + 5 * T**3. 

2 main => 

3 N = 11, 

4 As = to_array([read_real() : I in 1..N]), 

5 foreach (I in N..-1..1) 

e Y = f (As [I] ) , 

7 if Y > 400 then 

s printf ("7 0 w TOO LARGE\n", I - 1) 

9 else 

10 printf ("7 0 w %w\n", I - 1, Y) 

11 end 

12 end. 

Listing 1.1. TPK algorithm in Picat 

Line 1 defines a function to calculate the value of / (a function in Picat is 
a special kind of a predicate that always succeeds with a return value). Lines 
2- 12 define the main predicate. Line 4 uses list comprehension to read 11 space- 
separated real numbers into array As. Line 5 defines a header of foreach loop: 
I goes from 11 to 1 with the step -1 (in Picat array indices are 1-based). Lines 
6-11 calculate the value of y and print the result using an ‘if-then-else’ construct, 
printf is similar to the corresponding C language function; °/ 0 w can be seen as a 
“wildcard” control sequence to output values of different types. 


2 The Problems 

For this section we have chosen a set of GCJ problems from different years to 
demonstrate different useful aspects of Picat: constraint programming, top-down 
dynamic programming with tabling, and the planner module. 


Triangle Area^l 

“Triangle Areas” is a problem from the round 2 of GCJ 2008. The problem gives 
integers N, M and A and asks to find any triangle with vertices in integer points 
with coordinates 0 < Xj < N and 0 < yi < M that has an area of ^, or to decide 
that it does not exist. 

“Triangle Areas” is almost perfect for solving with constraint logic program¬ 
ming. Variables are discrete, constraints are non-linear, and we are looking for 
any feasible solution. To come up with an effective model we need to notice 
that one vertex of the triangle can be chosen arbitrarily. With this observation, 
the most convenient way to calculate the doubled triangle area is to place one 
vertex in (0, 0); then 2 S = \x 2 y 3 — X 3 IJ 2 1. (The same formula can be used in an 
imperative solution.) 

For this problem we present complete source code of the solution. For sub¬ 
sequent problems we omit the main predicate to save space. 

5 Problem link: http://goo.gl/enHWlq 



1 import cp. 

2 import util. 

3 model(N, M, A, Points) => 

4 [X2, X3] :: 0..N, 

5 [Y2, Y3] :: 0..M, 

e A #= abs(X2 * Y3 - X3 * Y2), 

7 Points = [X2, Y2, X3, Y3]. 

8 do_case(Case_num, N, M, A) => 

9 printfC'Case #°/,w: ", Case_num) , 

10 if model(N, M, A, Points), solve(Points) then 

11 printfC'O 0 °/ 0 s\n" , join( [to_string(V) : V in Points])) 

12 else 

is printlnC"IMPOSSIBLE") 

14 end. 

is main => 

16 C = read_int(), 

17 foreach (Case_num in 1..C) 

is N = read_int(), M = read_int(), A = read_int(), 

19 do_case(Case_num, N, M, A) 

20 end. 

Listing 1.2. Complete Picat program for the “Triangle Areas” problem 

Lines 1-2 load Picat modules for constraint programming and utility func¬ 
tions. Lines 3-7 define the model with input parameters N, M, A and a list of 
output parameters [X2, Y2, X3, Y3]. : : and #= are from the ‘cp’ module. 
With : : we define possible domains for X2, X3, Y2, Y3 variables, and #= from 
‘cp’ constraints both left and right parts to be equal. After model evaluation X2, 
X3, Y2, Y3 variables will not necessarily be instantiated to concrete values, but 
they will have reduced domains with possible delayed constraints and will be 
instantiated later with solve. 

Lines 8-14 define the do_case predicate to process a single input case. Line 9 
outputs case number according to the problem specification. Lines 10-14 are an 
‘if-then-else’ construct that outputs point coordinates if it is possible to satisfy 
our model predicate and solve (assign concrete values from the domain to every 
variable) the constraint satisfaction problem, or “IMPOSSIBLE” otherwise. Line 
11 uses an interesting Picat feature - list comprehension - which is very similar 
to what Python and many other modern programming languages have. 

Lines 15-20 define the main predicate that reads the number of test cases C 
and for each test case reads N, M, A parameters and executes do_case. 

This Picat program is very similar to our constraint programming solution 
in ECL’PS® CLP [Tj. 

A possible imperative solution in a mainstream programming language re¬ 
quires a more in-depth analysis of the problem. First observations will be the 
same. We will also note that it is impossible to find required triangle if A > 
M x N, and for A = M x N triangle (0, 0), ( N , 0), (0, M ) is a valid answer. Now, 
for A < MxN we can represent A as M(A div M) + (A mod M), 0 < A div M < 


N, 0 < A mod M < M. If we match this representation with the area formula, 
we can see that points (0, 0), (1, M), and {—A div M, A mod M) form a triangle 
with area j-- If we shift this triangle A div M units in positive direction along 
the x axis, we will get a triangle {A div M, 0), (A div M + 1, M), (0, A mod M) 
that will match all the requirements. 

Arguably, our declarative solution in Picat is simpler and leaves less space 
for a possible mistake. 

Interestingly, in our tests the running time of our solution in Picat 0.9 on 
small input is about 2.5 times larger than on the large input (table [l]). This is 
probably related to the implementation details and could change in the future 
versions. 


Welcome to Code Jair0 

“Welcome to Code Jam” is a problem from the qualification round of GCJ 2009. 
The task is to calculate the last 4 digits of the number of times the string 
“welcome to code jam” ( S ) appears as a subsequence of the given string (T). 

This is a typical dynamic programming problem. The problem state dp[i] [j] 
is the number of times the substring of T of length i contains the substring 
of S of length j (modulo 10000). The recurrence relation is: if T[i] = S'fj], 
dp[i\\j] = dp[i - 1 ]\j - 1] + dp[i - 1][j], otherwise dp[i\\j] = dp[i - l][j]. 

Our Picat program uses tabling [3! (a kind of memoization) to implement 
the described dynamic programming solution in a top-down fashion. 

1 s() = to_array("welcome to code jam"). 

2 table 

3 ways(_, _, 0) = 1. 

4 ways(_, 0, _) =0. 

5 ways(T, I, J) = W => 

e S = s(), 

7 if T [I] == S [J] then 

8 W = (ways(T, I - 1 , J) + ways(T, I - 1 , J-l)) mod 10000 

9 else 

10 W = ways(T, I - 1 , J) 

11 end. 

12 do_case(Case_num, T) => 

13 W = ways(T, length(T), length(s())), 

14 printf("Case #%w: % 04 d\n", Case_num, W). 

Listing 1.3. Picat solution for the “Welcome to Code Jam” problem 

Line 1 defines the string S from the problem statement as a functional fact. 
Lines 2-11 defines recursive ways function for dynamic programming. The calls 
to this function are automatically tabled (memoized) because of the table decla¬ 
ration. Lines 3 and 4 describe base cases for the recursion. Lines 5-11 specify the 


Problem link: http://goo.gl/qeLls4 



recurrence relation. The do_case predicate in lines 12-14 calls the ways function 
and prints the results according to the problem specification. 

An imperative solution can rely on the same recurrence relation, but might 
require more code to implement it either as a bottom-up dynamic programming 
or as a top-down recursion with memoization. 


Bribe the Prisoner^ 

Bribe the Prisoners was the hardest problem from the round 1C of GCJ 2009. 
In it we have an array of P prison cells, each cell is either empty or contains 
a prisoner. Every time a prisoner from one of the cells is released, all prisoners 
housed on either side of that cell until cell 1, cell P, or an empty cell get one 
coin each. Initially all cells contain prisoners. Given a list of indices of prisoners 
to be released, find the minimum total number of coins that will be spent if the 
prisoners will be released in an optimal order. 

This is an another dynamic programming problem. For each pair of cells 
A < B, dp[A][I3] is the best answer if prisoners occupy only cells from A to 
B, inclusive. If the first prisoner between A and B to be released is in cell X , 
(. B — A) coins are to be paid out immediately after his release, and then the 
smaller subproblems dp[A][X — 1] and dp[X + 1][-B] have to be solved. The final 
answer dp[l][P] corresponds to the initial state of all cells occupied. 

Our Picat program uses mode-directed tabling [B;. 

1 table (+, +, +, min) 

2 cost(A, B, FreeList, Cost) ?=> 

3 foreach(X in FreeList) 

4 (X < A ; X > B) 

5 end , 

e Cost = 0 . 

7 cost(A, B, FreeList, Cost) ?=> 

8 member(X, FreeList), 

9 X >= A, X =< B, 

io cost(A, X - 1 , FreeList, CostLeft), 

n cost(X + 1 , B, FreeList, CostRight), 

12 Cost = B - A + CostLeft + CostRight. 

13 do_case(Case_num, P, FreeList) => 

14 cost(l, P, FreeList, Cost), 

is printfC'Case #°/ 0 w: 7 „w\n", Case_num, Cost). 

Listing 1.4. Picat solution for the “Bribe the Prisoners” problem 
The first line declares the tabling mode for the cost predicate: first 3 pa¬ 
rameters are input parameters, and the last parameter is an output parameter 
(the cost of releasing all prisoners in FreeList that occupy cells in the [A; B] 
range) that must be minimized. Lines 2-12 define two clauses of the cost pred¬ 
icate; both clauses are non-deterministic, and this is stated by using ?=> syntax 

7 Problem link: http://goo.gl/pSbrTk 



instead of =>. The first clause states that if no prisoners in FreeList occupy 
cells between A and B, the cost of their release will be 0. The second clause 
calculates the release cost of prisoner X as described by the recurrence relation. 
The do_case predicate in lines 13-15 calls the cost function for the whole range 
of cells and prints the result according to the problem specification. 

As with the previous problem, an imperative solution can use the same re¬ 
currence relation, but might require more code for a bottom-up or top-down 
approach implementation, including explicit comparison of release costs of dif¬ 
ferent prisoners to find the minimum. Our Picat solution replaces most of the 
auxiliary code with a single table declaration. 


Osmo^l 

“Osmos” is a problem from the round IB of GCJ 2013. The problem describes 
“motes” of different integer sizes. One mote (Armin) is controlled by a player, 
the rest are passive. If Armin is of size A', it can absorb any passive mote of size 
Y < X and grow to size X + Y as a result. You are given the initial size of Armin 
and the sizes of passive motes. You can add a passive mote of any positive size, or 
you can remove any existing passive mote. Minimize the number of addition and 
removal operations required for Armin to be able to absorb all passive motes. 
Out Picat program uses the planner module [7]. 

To come up with an effective planning solution we need to notice that there 
always exists an optimal solution in which Armin absorbs passive motes in order 
from smallest to largest (if there is a pair of motes absorbed in different order, 
they can be swapped without increasing the number of operations needed). 

1 import planner. 

2 final([_, []]) => true. 

3 action![Armin, Others], NewState, Action, Cost) ?=> 

4 Others = [Min I Rest], 

5 Armin > Min, 

6 NewArmin is Armin + Min, 

7 Action = absorb, 

8 Cost = 0 , 

9 NewState = [NewArmin, Rest]. 

10 action([Armin, Others], NewState, Action, Cost) ?=> 

11 Others = [Min I _Rest], 

12 Armin =< Min, 

13 append(NewOthers, [_], Others), 

14 Action = remove, 

is Cost = 1 , 

1 6 NewState = [Armin, NewOthers]. 

17 action![Armin, Others], NewState, Action, Cost) ?=> 

is Others = [Min I _Rest] , 


Problem link: http://goo.gl/0N5zB8 



19 Armin =< Min, 

20 Newltem is Armin - 1 , 

21 NewOthers = [Newltem I Others], 

22 Action = add, 

23 Cost = 1 , 

24 NewState = [Armin, NewOthers]. 

25 do_case(Case_num, Armin, Others) => 

26 Limit = length(Others), 

27 best_plan([Armin, sort(Others)], Limit, _Plan, Cost), 

28 printfC'Case #%w: °/ 0 w\n" , Case_num, Cost). 

Listing 1.5. Picat solution for the “Osmos” problem 

Solving a planning problem in Picat requires a final predicate and an action 
predicate. Line 2 defines the final predicate which has one parameter - the 
current state - and succeeds if the state is final. In our program a state is 
represented as a 2-element list: the first item is the Armin size, and the second 
item is a sorted list of the sizes of passive motes (Others). A state is final if the 
Others list is empty. 

Lines 3-24 define the action predicate which has three clauses - one for 
absorb, remove and add actions - and has four parameters: current state, new 
state, action name, and action cost. Lines 3-9 define the absorb action which 
can be used if Armin is bigger than the first of the other motes at the cost of 0. 
Lines 10-16 define the remove action which removes the last (the largest) mote 
from Others at the cost of 1. The append predicate and the [ I ] syntax for 
getting the head and the tail of a list work exactly the same way as in Prolog. 
Lines 17-24 define the add action which adds a mote of size Armin — 1 to the 
beginning of the Others (so it can be absorbed by the next absorb action) at 
the cost of 1. 

Picat’s predicate for finding an optimal plan best_plan has two input 
parameters: the initial state and the resource limit, and two output parameters: 
the best plan and its cost. To find an optimal plan the system uses tabling and 
iterative deepening depth-first search-like algorithm. If no plan was found and 
the maximum resource limit was reached, the predicate fails. In this problem 
the resource limit for best_plan is the initial number of passive motes, because 
there is an obvious plan of this cost to remove all the motes. 

This solution is a declarative specification of the problem statement which 
relies on just a few observations about the problem. An imperative solution 
would require much more insight into the problem. One could notice that in an 
optimal solution if a mote is removed, all motes of equal or greater sizes are also 
removed (if one of larger motes is absorbed, so can be the mote in question). 
Thus, a greedy solution is: keep absorbing passive motes from smallest to largest 
while absorbing the next one is possible. After this, either remove all passive 
motes left or keep adding motes of size one less than Arinin’s current size and 
immediately absorbing them until Armin can absorb the smallest passive mote 
left. Repeat until Armin absorbs the last of the given passive motes. 



Table 1. Running times for small (4 minutes time limit) and large (8 minutes) input® 


Problem 

Technique 

Small Large 

Triangle Areas 

constraint programming 

2.4s 

0.9s 

Welcome to Code Jam dynamic programming 

0.0s 

0.3s 

Bribe the Prisoners 

dynamic programming 

0.0s 

4.7s 

Osmos 

planning 

0.0s 

0.1s 


3 Conclusions 

We have given several examples of declarative solutions for GCJ problems with 
Picat using constraint logic programming and tabled logic programming. 

We considered using Picat’s mixed integer programming module which might 
be useful for solving many GCJ problems [T], but currently there is no easy way 
to suppress log messages written to standard output by the underlying solver. 

Compared to Prolog, Picat code can be more compact because of functions 
(function calls can be nested, so there is no need for intermediate variables), list 
comprehensions, and more convenient console input/output. Also, while many 
modern Prolog systems have loop constructs, Picat loop syntax looks much 
cleaner because neither global nor local variables need to be explicitly declared. 

Running times of our Picat programs are several orders of magnitude smaller 
than the time limit imposed by GCJ rules (table []}. 

We also have found that GCJ problems can be complex and large enough to 
exercise many different aspects of a programming language implementation: we 
discovered and reported two serious bugs in the version 0.8 of the Picat system 
while working on this paper (the bugs were promptly fixed for Picat 0.9). 
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Results were obtained on a 64-bit Linux machine with Intel Core i7-4900MQ CPU 
@ 2.80GHz and 16GB RAM using Picat 0.9. 






